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ABSTRACT 

We present models of the extragalactic background light (EBL) based on several 
scenarios of galaxy formation and evolution. We have treated galaxy formation with 
the Press-Schecter approximation for both cold dark matter (CDM) and cold+hot dark 
matter (CHDM) models. Galaxy evolution has been treated by considering a variety 
of stellar types, different initial mass functions and star formation histories, and with 
an accounting of dust absorption and emission. We find that the dominant factor 
infiuencing the EBL is the epoch of galaxy formation. A recently proposed method for 
observing the EBL utilizing the absorption of ~ 0.1 to 10 TeV gamma-rays from active 
galactic nuclei (AGN) is shown to be capable of discriminating between different galaxy 
formation epochs. Observations of TeV gamma-ray bursts (GRB) would constrain or 
eliminate models in which the GRB sources lie at cosmological distances. 
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1. Introduction 

We have developed an indirect test of the era of galaxy formation by modeling the extragalactic 
background light (EBL) in the near infrared portion of the spectrum — a domain difficult to 
observe directly due to galactic and zodiacal contamination (Matsuura et al. f994, Mattila f990). 
A new method for probing this spectrum has been suggested (Stecker et al. f992) which involves 
the absorption of TeV 7-rays by the EBL. TeV 7-rays have recently been detected from Markarian 
421 (Punch et al. 1992), an AGN at a redshift of 2: = 0.031. If several more distant TeV sources 
exist, it should be possible to observe for the first time the EBL, and thus constrain the timescale 
over which galaxies form. The potential of this technique is illustrated by the constraints it places 
upon an earlier claim for a detection of the EBL (Matsumoto et al. 1988). 

It has been known for many years that high energy 7-rays from sources at cosmological 
distances will be absorbed by a diffuse background of long wavelength photons (Gould & Schreder 
1967) through electron-positron pair production. This point had been of little interest, however, 
until the discovery of a class of high energy 7-ray emitters by the EGRET detector aboard 
the Compton Gamma Ray Observatory (Lin et al. 1992). These sources, all AGNs of various 
subclasses, were found to emit a substantial amount of energy in the GeV energy domain with 
a roughly E~'^ spectrum. While the physical description of the emission process is uncertain, it 
is possible that this spectrum could continue into the TeV range for some of these sources, a 
hypothesis given support by the detection of Mrk 421 at TeV energies. Mrk 421 is the nearest 
of the EGRET sources and should suffer from only limited EBL absorption of its 7-rays as we 
will show, but sources at greater distances should show evidence of absorption in the TeV domain 
which can serve as an indicator of the era of galaxy formtion. 

2. Extragalactic Background Models 

Models of the EBL require several inputs: a library of stellar spectra and a function describing 
the number distribution by mass of a stellar population (the initial mass function, IMF); a 
characterization of galactic classes (elliptical, spiral) and the evolution of their respective stellar 
populations (the star formation rate, SFR) and dust content; a function describing the formation 
of galaxies as a function of redshift; and assumptions about the geometry of the Universe — the 
density parameter (rio)j the cosmological constant (A) and the Hubble constant {Hq). Such 
models have been created before (Partridge & Peebles 1967; Yoshii & Takahara 1988; Franceschini 
et al. 1994) with a variety of different assumptions, but each has a perhaps oversimplified model of 
galaxy formation. The conclusion reached by these earlier modelers is that the EBL proves to be 
a poor test of cosmology, meaning geometry, as the uncertainties in the IMF and SFR overshadow 
differences due to and Hq. We concur that the EBL proves a poor test of geometry. But this 
by no means limits its usefulness as a probe of cosmology, for the dominant factor infiuencing the 
EBL is the era of galaxy formation, which in modern theories depends strongly on the nature of 
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dark matter. For simplicity, we have chosen to consider here only those models with a flat {k = 0) 
geometry and A = 0, i.e. Q = 1; the arguments that the age of the Universe to ^ 13 Gyr then 
require Hq ~ 50 km s~^ Mpc~^. Considering also < 1, especially with A > 0, would tend to 
increase differences due to the galaxy formation epoch. 

Observational work has suggested several models for the IMF; we have considered the full 
range predicted by various authors (Scalo 1986; Basu & Rana 1992). The SFR has been given the 
functional form of a decaying exponential in time (Bruzual & Chariot 1993). The instantaneous 
recycling approximation has been used to model the evolution of the gas content and metallicity 
within each galaxy (Tinsley 1980). Dust absorption has been treated in a fashion similar to 
Guiderdoni & Rocca-Volmerange (1987); the dust emission spectrum is taken from Desert et 
al. (1990). Galaxy formation has been treated with the Press-Schecter approximation, which 
expresses the number density of galaxy halos of a given mass as a function of redshift based 
on the results of N-body simulations (Klypin et al. 1993). Two different cosmological models 
for structure formation have been considered: the cold dark matter (CDM) model, representing 
a moderately early era of galaxy formation, and cold+hot dark matter (CHDM), representing 
a late era of galaxy formation. The latter is normalized to COBE and both are based on the 
simulations of Klypin et al. (1993). There is a difficultly in describing a galactic population in this 
manner in that the Press-Schecter function merely predicts the abundance of dark matter halos 
and makes no distinction between a cluster of galaxies and a single galaxy. One must therefore 
also assume a maximum galactic mass; we set the upper limit at M^ax = 5 X 10^^ Mq based on 
a comparison between the galactic luminosities obtained for a power law mass to light ratio and 
the normalization process described below, and we have checked that the results do not depend 
sensitively on M^ax- 

The flnal parameter to be set is the normalization of the SFR, or the efficiency of star 
formation for a given galactic type. This can be directly tied to the local galaxy population 
through a comparison of model galactic luminosities with those of the nearby population. A 
number of studies have been made of the galaxy local luminosity function (LLF). Most recent is 
the work of Loveday et al. (1992) based on the Stromlo-APM southern sky survey (LEPM LLF), 
the LLF of Efstathiou et al. (1988) based on an average of several surveys (EEP LLF), and the 
LLF of DeLapparent et al. (1989) based on the CfA survey (LGH LLF). Each has characterized 
the luminosity distributions of the galaxies in the blue magnitude band as a Schecter function, 

(f){y) dy = (f)^y"e~y dy where y = L/L^b- (1) 

If one makes the assumption that the most massive galaxies are also the most luminous, one can 
compare the number densities of the LLF and the Press-Schecter function (at z = 0) and require 
that the average age model galaxies, at the present epoch, have the expected luminosities by 
deriving a function which expresses the SFR normalization as a function of galactic mass. The 
luminosity of the galaxies is thus normalized by forcing the model galaxies to duplicate the local 
population. 
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The model results are shown in Figure 1. Note the large separation between the galaxy 
formation model predictions; this will be seen for any combination of model parameters given that 
we include both early and late galaxy formation models and that we require that the present day 
model galaxies match the local population. The larger EBL flux predicted by the CDM model 
arises for several reasons: (1) stars have been contributing their light to the EBL for ~ 2 Gyr 
longer than in the CHDM case, (2) the initial burst of star formation in early-type galaxies has 
redshifted from the optical to the near IR, and (3) the galaxies are older at a given redshift and 
hence are composed of more evolved stars, producing a brighter flux in the red and near IR. In 
addition to the parameter choices discussed above, we have also varied the other quantities not 
flxed by observational data to ensure that no large deviations were found. (A more detailed paper 
describing our models and results is in preparation.) 



The cross section for the absorption of 7-rays through e"'"e pair production is given by Gould 
& Schreder (1967). This cross section is maximized when 



where e is the EBL photon energy and E is the 7-ray energy. The photon number density for 
the EBL models drops off above 1 eV (~ 1 //m) and below 0.1 eV; the absorption effect due to 
the stellar component of the EBL will be most pronounced for 7-ray energies of 300 GeV to 3 
TeV — see Figure 2. Above 10 TeV, it is likely that the dust-emitted component of the EBL makes 
the Universe mostly opaque, and above 100 TeV the cosmic microwave background completely 
absorbs any 7-rays from sources at cosmological distances. There is also an effect due to the 
redshift of the source. The absorption at some redshift along the path from the source will occur 
between a background photon of energy eo(l -|- z) and a 7-ray of energy £'0(1 + z), hence the 
lower cutoff in a source's spectrum will scale with ~ {1 -\- z)~'^. For sources at z > 0.5 there is an 
additional effect from the galaxy formation models in that the magnitude of the EBL at large z is 
quite small for cosmologies with late galaxy formation such as CHDM. 

This last point, combined with the exponential dependence of 7-ray absorption upon the 
model predictions of the EBL, gives this test its potential power — provided that there are several 
sources which do possess a E~'^ 7-1'ay spectrum stretching into the TeV range, and that there are 
instruments with the sensitivity to detect fluxes on the order of a few events per cm^ day. The 
one source known to emit TeV 7-rays, Mrk 421, may already constrain an earlier measurement of 
a diffuse infrared emission component of the sky which was attributed to the EBL — see Figure 3. 
The Mrk 421 spectrum also shows some evidence for a cutoff above 2 TeV; the magnitude of this 
cutoff is too large to explain by absorption with physically motivated EBL models and may instead 
represent absorption intrinsic to the source. This explanation would suggest that the mechanism 
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responsible for the production of the 7-rays was based on an electron Comptonization process 
(Dermer & Schlickeiser 1993; Sikora et al. 1994) rather than a proton based model (Mannheim 
1993). 

The EGRET team has observed at least 31 AGNs to date in the GeV energy range (Fitchel 
et al. 1994) of which 10 are at 2: < 1 and all with fluxes greater than that of Mrk 421. As of yet, 
only Mrk 421 has been seen in TeV light; the Whipple Observatory group has published upper 
limits on the TeV emission for two of the nearer GeV sources (3C 273 and 3C 279) (Fennell et 
al. 1993), but nothing conclusive can yet be said. EGRET has also observed several 7-ray bursts 
(GRB) in GeV light (Fitchel et al. 1994) with a comparable spectrum to the AGNs. It is possible 
that the process responsible for the GRBs may also produce a TeV flux, and this flux could be 
used to probe the extragalactic background as well. Given the strong absorption for < 200 
GeV shown in Figure 2 for all of our models for sources at z > 0.5, the identiflcation of any GRB 
at TeV energies would serve to constrain, or even eliminate, most cosmological theories of the 
mechanism responsible. 

Observations at 100 GeV, below the cutoff expected for any model of the EBL, should 
help clarify whether the spectrum of a given source shows intrinsic absorption. The Whipple 
telescope is currently sensitive down to ~ 500 GeV. There are a number of instruments planned 
(and some funded) which can perhaps achieve lower thresholds as well as several other detectors 
to complement Whipple, chiefly the GLAST satellite, sensitive from 10 MeV up to ~100 GeV 
(GLAST workshop, 1994); the MILAGRO detector, operating above 300 GeV; and the Themis and 
Solar One arrays, operating at energies ^^100 GeV (Rene Ong 1994, personal communication). If 
further sources are found to emit TeV 7-rays, the new detectors, combined with existing facilities, 
should be able to constrain many theories of galaxy formation. There are also a number of 
detectors operating at energies above 10 TeV— the CASA-MIA, CYGNUS, HEGRA, and TIBET 
arrays (Cronin et al. 1993); given the absorption of such 7-rays predicted in this work, it unlikely 
that such detectors will see any 7-rays from sources at cosmological distances. 

We have beneflted from encouragement and helpful conversations with David Dorfan, Spencer 
Klein, Richard Lamb, Rene Ong, Floyd Stecker, Trevor Weekes and especially Anatoly Klypin. 
This research was supported by NSF REU funds and a grant from the DOE at UCSC. 
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Fig. 1. — The EBL at z = for a variety of models. The upper set of curves is for the CDM model 
of galaxy formation, the lower set for CHDM. The CDM model has ^CDM = 0.9, = 0.1 and 
linear bias of 1.5; the bulk of the galaxy formation takes place at 2: = 1 to 3. The CHDM model 
is based on ^CDM = 0.6, fij, = 0.3 and Qb = 0.1; the bulk of galaxy formation takes place at a 
z = 0.2 to 1. Star formation occurs at a rate cx e~*/^ where r = 0.5 Gyr for elliptical galaxies 
(28% of the total) and r = 6 Gyr for spirals (72%). The number distribution of stars by mass 
was taken to be a power law of the form N cx m^~^; three values for F were considered for each 
galaxy formation model. Two of those presented here represent the extreme values for both CDM 
and CHDM galaxy formation models, with Model 1 representing the bluest spectra and Model 3 
the reddest; Model 2 is midway between these extremes. The dust models are based on a three 
component model corresponding to the PAH molecules (~ 10 to 30 /um), warm dust from active 
star forming regions (30 to 70 /um), and cold "cirrus" dust (70 to 1000 ^um). The near vertical line 
at the right represents the CMBR. Data of Yoshii & Takahara (1988) is taken from their model with 
cosmological parameters similar to those used herein; the Tyson (1990) data is a direct estimate of 
the EBF based solely on number counts of galaxies. The DeJager et al. (1994) model is based on 
a power law EBF fit to the possible absorption of Mrk 421 7-rays shown in Fig. 3. Observations 
(triangles) and upper limits (arrows) are from Matsumoto et al. (1988); a similar experiment failed 
to confirm this measurement (Matsuura et al. 1994). 
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Fig. 2. — The expected 7-ray flux from a source at three different redshifts {z = 0.03, 0.2, 0.5) with 
a spectrum given by dN/dE = 5 X IO'^Eq'^y cm~^ s~^ GeV~^ (typical of the EGRET detected 
AGNs). The absorption is calculated by assuming an exponential attenuation where the optical 
depth is given by 

t{E) = — / / -Vl + zn{e, z)a{xEe{l + z)"^) dz dx de 

Ho Jo Jo Jlm^c^ /Ex(l+z)^ 2 

where z^ is the redshift of the source; x = 1 — cos 9, where 9 is the encounter angle of the photons; 
n{e,z) is given by the EBL models in Fig. 1; and a is the pair production cross section. For each 
redshift, the three lower curves within each set represent Models 1, 2 and 3 from top to bottom. 
The upper curve in each set shows the absorption due to the stellar component of the EBL for the 
Model 2 parameters. 
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Fig. 3. — The 7-ray flux from Mrk 421. Data is from the EGRET detector (Lin et al. 1992) (circles) 
and Whipple Observatory (Mohanty et al. 1993) (triangles). The Whipple data is somewhat 
uncertain due to difficulties with their energy calibrations (Richard Lamb and Trevor Weekes 1994, 
personal communications). The unabsorbed spectrum was determined by fltting a power law to 
the two data sets without assuming any absorption in the TeV domain; including absorption would 
lower the power somewhat (the EGRET data alone suggests an E~^'^^ spectrum) and bring both 
Model 2 curves into agreement with the flrst four Whipple points. The £-2.55 fj-Qj-,^ ^j^g Y)q Jager 
et al. (1994) model in Fig. 1. The 1 to 3 //m curve shows the absorption which would be expected if 
one interpets the measurement of Matsumoto et al. (1988) as a background of extragalactic origin; 
due to the sharp cutoff below 1 TeV, reasonable flts of the observed spectrum to this background 
are not possible. 



